Discovery of a Radio Source following the 27 
December 2004 Giant Flare from SGR 1806-20 
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Over a decade ago it was established that the remarkable high energy tran- 
sients, known as soft gamma-ray repeaters (SGRs), are a Galactic popula- 
tion 1,2 and originate from neutron stars with intense 10 15 G) magnetic fields 
("magnetars" 3 ). On 27 December 2004 a giant flare 4 (fluence 5 ^ 0.3 erg cm 2 ) 
was detected from SGR 1806—20. Here we report the discovery of a fading 
radio counterpart. We began a monitoring program from 0.2 GHz to 250 GHz 
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and obtained a high resolution 21-cm radio spectrum which traces the inter- 



vening interstellar neutral Hydrogen clouds. Analysis of the spectrum yields 
the first direct distance measurement of SGR 1806—20. The source is located 
at a distance greater than 6.4 kpc and we argue that it is nearer than 9.8 kpc. 
If true, our distance estimate lowers the total energy of the explosion and re- 
laxes the demands on theoretical models. The energetics and the rapid decay 
of the radio source are not compatible with the afterglow model that is usually 
invoked for gamma-ray bursts. Instead we suggest that the rapidly decaying 
radio emission arises from the debris ejected during the explosion. 



On 3 January 2005 we observed SGR 1806— 20 with the Very Large Array 
(VLA) and identified and promptly reported 6 a new radio source at right ascension 
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«j2ooo =18 h 08 m 39.34 s and declination 5j 20 oo =-20°24'39.7 // (with an uncertainty of ±0.1" 
in each coordinate) coincident with the quiescent X-ray counterpart. 7 In Table 1 we report 
the results of a subsequent monitoring program undertaken with the VLA, the Giant Me- 
trewave Radio Telescope (GMRT), the Australia Telescope Compact Array (ATCA), the 
Nobeyama Millimeter Array (NMA) and the Institut de Radioastronomie Millimetrique 
(IRAM) 30m Telescope. 

The radio source decays in all frequency bands, but the behaviour is complex (Fig- 
ure 1). At each band we model the flux by a power law, S„(t) oc t a , but allow for changes 
in the temporal indices a ( "breaks" ) at two epochs. These breaks are clearly seen in our 
highest signal-to-noise ratio data. Following the first break (9 days, postburst) the light 
curve steepens to — 4 ^ a ^ —3. The radio source 8 from SGR 1900+14 following the 27 
August 1998 giant flare 9 showed a similar rapid decay at 8 GHz. Subsequently around 
day 14 the light curve flattens to a ~ — 1. At any given epoch, the radio spectrum can be 
modeled by a power law, S u oc iA The spectral index, ft, steepens with time, changing 
from ~ —0.7 to —0.9 (see Figure 1 and Supplemental Information). 

We confirm claims that the source is resolved 10 by an independent analysis. We find 
that it is elongated with a major axis 6 ~ 77 milliarcsecond (mas) and an axial ratio of 
2:1 (Table 2). We considered four expansion models, 6 cx t s with unconstrained s and 
plausible expansion models (s = 0, 2/5, and 1). The best fit model corresponds to no 
expansion (s = 0.04 ± 0.15). However, due to the limited range of our observations we 
prefer not to model the dynamics of the explosion. 

We took advantage of the brightness of the radio source and obtained a high resolution 
spectrum (Figure 2, second panel) centered around the 21-cm line of atomic Hydrogen 
(HI). Intervening interstellar clouds appear as absorption features in the spectrum. These 
clouds are expected to participate in the rotation of the Galaxy and the absorption features 
allow us to infer "kinematic" distance estimates. Such estimates have several complica- 
tions. First, in the inner Galaxy the radial velocity curve is double-valued (see Figure 2, 
third panel) leading to a "near" distance estimate (dj) and a "far" distance estimate (d u ) 
for each velocity. Second, in some directions, there are features with non-circular mo- 
tion e.g. the "3-kpc expanding arm" and the "—30 kms -1 spiral arm". 11 Finally, in the 
innermost part of the Galaxy there is a deficit of cold gas. 12 
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Significant HI absorption toward SGR 1806— 20 is seen over the velocity range —20 
kms" 1 to +85 kms -1 (Figure 2, second panel). There is also a weak (2.5-er) absorption 
feature coincident in velocity with a clearly detected 12 CO(1-0) emission feature identi- 
fied 13 as MC94 (Figure 2, first panel). Adopting a simple Galactic rotation curve with 
a circular velocity = 220 kms -1 and a Galactic center distance R =8.5kpc, the near 
distance to SGR 1806-20 (for V LSR =95 kms" 1 ) is d z =6.4kpc. 

The two HI emission clouds seen at velocities above 100 kms -1 toward SGR 1806— 20 
(Figure 2, first panel), with no corresponding HI absorption, may be used to infer an upper 
limit to the distance provided that we can be reasonably certain that cold neutral gas 
exists at these velocities. The HI absorption spectrum toward the nearby (A9 = 0.77°) 
extragalactic source J1811— 2055 shows a strong and broad absorption feature between 110 
and 130 kms -1 (Figure 2, fourth panel). The only HI emission in this direction 14 above 
60 kms -1 corresponds to an HI cloud at this same velocity. This feature can be traced in 
absorption toward several other extragalactic radio sources in this direction, 15 suggesting 
that cold gas at ~120 kms -1 is widespread. Adopting the same Galactic rotation curve 
as above, the absorbing cloud at +120 kms -1 can either be located at 7kpc or 9.8 kpc 
(see Figure 2, third panel). We thus suggest an upper limit to the distance, d n =9.8kpc. 

Our new distance estimate is smaller than previous (indirect) values 16 ' 11 of 12 to 15 kpc. 
Accepting our estimate has several important implications. It results in a reduction of the 
total energy released (oc d 2 ) as well as the rate of such events in nearby galaxies 17 (oc d 3 ), 
and calls into question the association 18 of SGR 1806— 20 with a star cluster along the 
same line-of-sight. Therefore, claims that magnetars originate from more massive stars 
than normal neutron stars 19 may be called into question. 

Next we consider the energetics of the material giving rise to the radio emission. As 
with many other radio sources, the power law spectrum can be attributed to energetic 
electrons with a power law energy distribution (dN/d'-f oc 7~ p ; here 7 is the Lorentz 
factor of electrons and we measure p = 2.24 ± 0.04 on day 7, which as is a typical value 
for strong shocks) which gyrate in a magnetic field and emit synchrotron radiation. We 
apply the minimum energy formulation for synchrotron sources 20 ' 21 to the radio spectrum 
(from 0.2 GHz to 100 GHz) of 3 January 2005 and find the energy of the radio emitting 
source and the associated magnetic field strength are U m in ~ lO 43 ^^^ 7 ^^ 7 erg and B min ~ 
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lU^ l7 e-^ mG; here, the distance is 10rf 10 kpc and the angnlar diameter is 75£ 75 mas 
(Table 2). 

Evidently the amount of energy released in the 7-ray flare, E Jjiso ^ 4 x 10 45 ergs _1 
(assuming unbeamed, isotropic emission), substantially exceeds U min . In contrast, the 
ratio C/ m i n /£' 7) i so is unity for GRBs and as a result the lower energy and longer duration 
emission is correctly regarded as arising from the shock of the circumburst medium (the 
"afterglow" model). Thus, based solely on energetics, there is no prima facie reason to 
suggest that the radio source is the afterglow of the 7-ray flare. 

Furthermore, as discussed above, the radio emission decays quite rapidly 9 days after 
the burst. Such a rapid decay is incompatible with the afterglow model (in the non- 
relativistic limit) for which we expect 22 a = 3(3 + 0.6. We conclude (in contrast to 
refs. 23, 24 and 25) that the radio emission must be powered by something very different 
from that which produced the 7-ray emission. 

In summary, the radio emission can be described by two components: a rapidly decay- 
ing component and a slowly decaying component. The latter becomes detectable when 
the former has decayed significantly. The rapid decay is phenomenologically similar to 
that seen from accreting Galactic sources (e.g. ref. 26) - the so-called "plasmon" model 
framework in which the radio emission arises from a ball of electrons and magnetic field 
which are initially shocked and then cool down by expansion. We make the specific sug- 
gestion that the radio emission up until about 2 weeks is a result of the shocking of the 
debris given off in the explosion (the "reverse shock"). In this framework the slowly de- 
caying component is the emission arising from the forward shock as the ejecta slams into 
the circumburst medium. A requirement of this suggestion is that the energy inferred in 
the slowly decaying component should be comparable to U ni \ n . Separately, we note that 
the comparable ratios U m i n / E ltiso and the temporal and spectral similarities of the giant 
flares from SGR1806— 20 and SGR 1900+14 suggest a common mechanism for launching 
these flares and similar circumstellar environments. 

Regardless of the suggestions and speculations, it is clear that radio afterglow is telling 
us something entirely different from that revealed by the 7-ray emission. If our suggestion 
of a reverse shock origin is correct then radio observations allow us to probe the ejecta. 
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Taken together it appears that rapid and intense radio monitoring of such flares will be 
highly fruitful in the future. 
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FIGURES & CAPTIONS 



Table 1: Flux density measurements of the transient radio counterpart to SGR1806— 20 from 
the VLA, GMRT, NMA, and ATCA as a function of frequency and time. The reported errors are 
1-0- . In addition to these measurements, we obtained IRAM-30m observations on 8 and 9 Jan- 
uary 2005 using MAMBO-2 at 250 GHz which show no detection with a value of 0.57 ±0.46 mJy 
at the position of the radio source. Finally, we detect linearly polarized emission from the source 
at the 1.5% to 2.5% level. See the Supplemental Information for observational details. 
a The epoch of the flare, t , was 27.90 December 2004. 
b ATCA observations in this column have a frequency of 8.6 GHz. 
c These values represent 2a upper limits. 

d The frequency is 1.06 GHz for the 16.37 January 2005 and 4.01 February 2005 GMRT obser- 
vations. 

Figure 1: Broadband temporal behavior of the transient radio source coincident with 
SGR1806— 20. The abscissa indicates days elapsed since the giant flare on 27.90 December 
2004. The displayed flux density measurements (denoted with symbols) were obtained in six 
frequency bands with the VLA, GMRT, and ATCA (Table 1). The error bars denote 1-cr uncer- 
tainties. With the exception of the 6.1GHz data (which is insufficiently sampled at early and 
late times and is not shown) , the light curves with v > 1 GHz are best fit by power-law models 
(shown as lines, S v oc t ai ) with two breaks at t\ ~ 9 days and ~ 15 days (see the Supplemental 
Information for exact values). The temporal index varies chromatically in the time before and 
after the first break (denoted by regions A and B respectively). The exponent value ranges from 
—2 <^ a a — 1 and —4 <J as <J — 3; here the subscript identifies the region of interest. After 
day ~ 15 (region C) the source decay flattens to ac ~ —0.9 at these frequencies which persists 
until day 51. Region B, the period of steep light curve decline, is shaded gray. The light curves 
with v < 1 GHz do not show these temporal breaks or late time flattening Apparently a single 
power law decay model with a = -1.57 ± 0.11 (0.24 GHz) and a = -1.80 ± 0.08 (0.61 GHz) 
provides a good statistical description of the data. 

Our substantial frequency coverage (over three decades) allows an excellent characterization 
of the spectrum. The spectrum is consistent with a single power law slope (S u cx v") at all 
epochs. On day 7, before the first temporal break, we find (3 = —0.62 ± 0.02. The spectrum 
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steepens to a value of f3 = —0.76 ± 0.05 (day 15), reaching (5 = —0.9 ±0.1 (days 21-51). 



Table 2: Source size measurements and 95% confidence limits of the radio source as measured 
with the VLA at 8.46 GHz. The source is clearly resolved at all VLA epochs. The best constraints 
on the source size come from the observations which occurred closest to the transit of the source 
on Jan 3rd and Jan 6th. The result is a source with size 6 ~ 77 mas with an axial ratio of ~ 0.5 
and a position angle (PA) of 60 degrees (measured clockwise from the North). The flux centroid 
did not change position within the limits of our astrometric accuracy (± 100 mas). 

The best fit model is consistent with no expansion, s = 0.04+0.15 (9{t) oc t s ) with a x 2 = 5.3 
with 4 degrees of freedom. Sedov- Taylor (s = 2/5) and free expansion (s = 1) model were also 
fit, and yield x 2 = 20.6 and x 2 = 91, respectively. These fits have five degrees of freedom. See 
the Supplemental Information for the details of the source size measurements. 

Figure 2: Cold atomic and molecular hydrogen spectra toward SGR 1806— 20. These spectra 
were used to derive a distance estimate for SGR 1806— 20. (Top Panel). HI emission (upper 
curve, thick line) in the direction of SGR 1806— 20, determined by averaging two adjacent spec- 
tra taken by Hartmann & Burton 14 at 1, 6=(10.0°,0.0°) and I, 6=(10.0° -0.5°). The lower curve 
(thin line) is the 12 CO(l— 0) spectrum (from ref. 11). For display purposes the brightness tem- 
perature has been scaled up by a factor of 11.4. (Second Panel). The HI absorption spectrum 
taken toward SGR 1806— 20. The two horizontal bars illustrate the radial velocity measure- 
ments 16 ' 11 of the nearby star LBV 1806-20 (36+lOkms- 1 and 10 ±20 kms" 1 ). The absorption 
spectra were made with the Very Large Array on 4 January 2005, using a 1.56 MHz bandwidth 
in both hands of polarization centered 50 kms -1 with respect to the local standard of rest. The 
bandwidth was divided into 256 channels each 6.1 kHz in width, or a velocity resolution of 1.3 
kms -1 covering a velocity range of —115 kms -1 to +215 kms -1 . (Third Panel). The distance 
as a function of radial velocity adopting a simple Galactic rotation curve with a circular velocity 
G = 220 kms -1 and a galactic center distance R =8.5kpc. (Bottom Panel). The HI absorption 
spectrum of nearby extragalactic source J1811— 2055 at /, 6=(9.8°,— 1.0°). 

The lower limit to the distance is firmly established by the 95 kms -1 absorption feature 
from MC94 (see Text). An upper limit to the distance of the SGR is suggested by the absence 
of strong absorption at +120 km s _1 , seen toward J1811— 2055 and several other extragalactic 
radio sources in this direction. 15 One could argue that the 120 km s -1 cold cloud is small and 
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not present along the line-of-sight to the SGR 1806—20. However, this hypothesis also requires 
the absence of any other cloud between DSkms" 1 (distance of 6.4 kpc) and 86kms _1 (distance 
of 10.6 kpc). The mean absorption coefficient drops in the inner Galaxy (R <^ 4.5) kpc, giving 
a mean free path between clouds of 2.3 kpc. 15 The distance interval from 6.4 kpc to 10.6 kpc 
corresponds to ~ 1.8 mean free paths. So, the probability of finding no clouds in this gap is 
16%. Thus our upper limit of 9.8 kpc is not a certainty but quite likely. 
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Epoch 


Beam 


Beam FA 


Fit Major Axis 


Axial Hatio 


Fit FA 




(mas) 


(degrees) 


(mas) 




(degrees) 


03 Jan 2005 


349 x 170 


11.8 


78.233 


34+ uyl 




04 Jan 2005 


593 x 173 


-40.3 


72.4+it 5 


0.00+°' 90 


691?? 


05 Jan 2005 

06 Jan 2005 

07 Jan 2005 
10 Jan 2005 


397 x 168 
329 x 178 
532 x 178 
560 x 161 


-25.7 
-16.5 
40.4 

-38.8 


<55+ 18 
oo -10 

-c 7+3.0 
'°-'-3.0 
78+26 
'8-18 
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u - DU -0.60 
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SUPPLEMENTAL INFORMATION 



Observational Details: For the VLA observations we used the standard continuum mode 
with 2 x 50 MHz bands, with the exception of the 1.46 GHz observation of 4 January 2005, which 
was taken in spectral line mode with 8 channels of width 3.1MHz. We used the extragalactic 
source 3C286 (J1331+305) for flux calibration, while the phase was monitored with J1820-254, 
J1751— 253, and J1811— 209. The listed flux densities and uncertainties were measured from the 
resulting maps by fitting an elliptical Gaussian model to the afterglow emission. The GMRT 
observations were performed in dual frequency mode with 16 MHz bandwidth divided into a 
total of 128 frequency channels for 610 MHz observations, and 6 MHz of bandwidth divided into 
64 channels for the 235 MHz observations. The observations at 1060 MHz at GMRT were carried 
out with a bandwidth of 32 MHz. 3C48 (J0137+331) and 3C 286 were used as flux calibrators, 
and J1822— 096 was used as the phase calibrator. These sources were also used for bandpass 
calibration. We obtained the flux densities of the source by fitting a Gaussian with a back- 
ground level plus a slope and removed the contribution from a nearby weak source. Because 
of the high density of radio sources in the galactic plane in which SGR 1806— 20 is located, the 
antenna system noise temperature has notable contributions from the sky within the telescope 
beam. This reduces the signal to noise ratio and an appropriate correction must be made to 
the observed flux (especially at low frequencies), since the flux calibrators which establish the 
flux scale lie well outside the galactic plane and are in an environment of less sky temperature. 
We applied a T sys correction factor for 3C48 of 3.88 and 1.93, and for 3C 286 of 3.87 and 1.8, 
for the 235 MHz and 610 MHz, respectively. Both the VLA and GMRT data were reduced and 
analyzed using the Astronomical Image Processing System. The ATCA observations were per- 
formed in snapshot mode with 100MHz of effective bandwidth. The amplitude calibrator was 
J1934-638, whereas J1711-251, J1817-254, and J1811-209 where used as phase calibrators. 
The last of these was observed in a rapid (3 minute) cycle mode to compensate for its poor 
phase stability. The flux densities were determined by performing a local parabolic fit to the 
peak closest to the known position of the source. The NMA observations were performed at 
102 GHz in D-configuration (the most compact configuration) on 4 January 2005, and in AB 
configuration (longest baseline configuration) on 12, 13 January 2005. We used NRAO530 for 



16 Cameron et al. 

the phase calibration, and assumed it to have a flux density of 2.3 Jy. 

Details of Source Size Measurements: The source sizes were measured by modeling 
the calibrated visibilities with the model-fitting procedure in DIFMAP. This procedure employs 
the Levenberg-Marquardt non-linear least squares minimization technique while fitting a 6 pa- 
rameter elliptical Gaussian to the visibilities. The errors were determined with DIFWRAP using 
the following scheme: the source size parameters were stepped in small increments around their 
best-fitted value to form a grid of values. At each grid point the source size parameters were 
held fixed while the other model parameters were allowed to 'relax' with 4 model-fitting rounds. 
The 95% confidence limits were determined by those models that had a Ax 2 < 12.8 as mea- 
sured from the best-fit total x 2 (Press, W. H., Teukolsky, S. A., Vetterling, W. T. and Flannery, 
B. P. Numerical Recipies in C. The art of scientific computing. Cambridge: University Press, 
2nd ed. 1992). As a check we used phase only self-calibration as well as phase and amplitude 
self-calibration, both of which give consistent source size measurements. We also used 30 second 
time-averaged data sets (to reduce the number of degrees of freedom), and found the best fit 
model parameters agreed to within the errors. 



Frequency 
(GHz) 


OLA 


ti 
(days) 


CtB 


*2 
(days) 


ac 


U.24U 


-1.7 ± U.I 










0.61 


-1.9 ±0.1 










1.4 


-2.0 ±0.2 


10.7 ±0.3 


-4.1 ±0.3 


13.8 ± 0.2 


-0.85 ± 0.2 


2.4 


-0.95 ±0.3 


9.8 ±0.2 


-3.5 ±0.2 


13.8 ±0.5 


-0.95 ± 0.2 


4.9 


-1.55 ±0.15 


8.8 ±0.2 


-3.1 ±0.2 


18.6 ±3.0 


-0.65 ± 0.3 


6.1 


-2.3 ±0.1 










8.5 


-2.00 ±0.15 


8.1 ± 0.3 


-2.8 ±0.24 


18.0 ±3.0 


-0.64 ±0.4 



Table : Summary of temporal indices and breaks in the light curve of SGR 1806—20 in 
7 frequency bands. The fits represent minimums in x 2 subject to the conditions that the two 
power-law slopes are continuous at the break point and disagree by more than 1-cr. The first 
and second break points are denoted by t\ and t2, respectively. The temporal decay indices 
(S u oc t ai ) are a a for t < ti, olb for ti < t < <2, and «c for t > ti- 



